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ABSTRACT 

Cosmic-ray (CR) electrons and nuclei interact with the Galactic interstellar gas 
and produce high-energy 7 rays. The 7-ray emission rate per hydrogen atom, called 
emissivity, provides a unique indirect probe of the CR flux. We present the measurement 
and the interpretation of the emissivity in the solar neighborhood for 7-ray energy 
from 50 MeV to 50 GeV. We analyzed a subset of 4 years of observations from the 
Large Area Telescope (LAT) aboard the Fermi Gamma-ray Space Telescope (Fermi) 
restricted to absolute latitudes 10° < |6| < 70°. From a fit to the LAT data including 
atomic, molecular and ionized hydrogen column density templates as well as a dust 
optical depth map we derived the emissivities, the molecular hydrogen to CO conversion 
factor Xco = (0.902 ± 0.007) x 10 20 cnL 2 (K km s -1 ) -1 and the dust-to-gas ratio 
X oust = (41.4 ± 0.3) x 10 20 cm' 2 mag -1 . Moreover we detected for the first time 
7-ray emission from ionized hydrogen. We compared the extracted emissivities to those 
calculated from 7-ray production cross-sections and to CR spectra measured in the 
heliosphere. We observed that the experimental emissivities are reproduced only if the 
solar modulation is accounted for. This provides a direct detection of solar modulation 
observed previously through the anticorrelation between CR fluxes and solar activity. 
Finally we fitted a parametrized spectral form to the heliospheric CR observations as 
well as to the Ferrm-LAT emissivity and obtained compatible local interstellar spectra 
for proton and Helium kinetic energy per nucleon between between 1 and 100 GeV and 
for electron-positrons between 0.1 and 100 GeV. 

Subject headings: cosmic rays - gamma rays: diffuse background - gamma rays: ISM - 
Galaxy: general - ISM: general - radiation mechanisms: non-thermal 


Introduction 


The theory behind cosmic rays and interstellar medium (ISM) interactions developed in the 
1950s (Hayakawa 1952 Hutchinson 1952 Ginzburg & Syrovatskij 1965 Stecker 1970[ 1971) was 
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first experimentally addressed with balloons and rockets experiments in the sixties. However the 
first evidence of correlation between high-energy 7 rays and Galactic structures was observed with 
the launch of the Third Orbiting Solar Observatory OSO-3 in March 1967 (jKraushaar et al.|l972). 


A 7 -ray H I emissivity of 1.6xl0 -25 ph atom -1 s -1 above 100 MeV was derived from comparison 
with an early H 1 21-cm line radiation radio survey. The Small Astronomy Satellite 2 (SAS-2) 
(Fichtel et al. 1975) launched in November 1972 collected 20 times more photons than OSO-3 and 


provided clearer evidence of correlation between high-energy 7 rays and hydrogen (Kniffen et al. 
1977). Lebrun & Paul (1979) compared the diffuse gamma-ray intensities measured by SAS-2 with 
atomic hydrogen (H 1 ) column density and derived an integrated H I 7 -ray emissivity above 100 MeV 
of (2.9±0.6)xl0 ~ 25 ph atom -1 s -1 . Later Lebrun et al. (1982) extracted the 7 -ray emissivity of H 1 


in three energy bands from the Cosmic Ray Satellite (COS-B) observations (Bennett et al. 1976). 


They assumed a power-law spectrum for the 7 rays and obtained a differential emissivity of 2.1 x 
10 -23 X E -2 ph atom -1 s -1 MeV -1 and an integrated flux above 100 MeV of 2.1± 0.3x 10 -25 ph 


atom 1 s 1 . Using the same data and a model for the interstellar emission that included 7 rays also 


originating from the molecular hydrogen and inverse-Compton scattering (IC), Strong et al. (1988) 


investigated the radial variation of the emissivity and the systematic uncertainties associated to the 


model. The Energetic Gamma Ray Experiment Telescope (EGRET) (Thompson et al. 1993) on 
the Compton Gamma-Ray Observatory, launched in 1991, provided a higher angular and energy 
resolution and collected about 4 times more 7 rays than COS-B above 100 MeV. In Strong & 


Mattox (1996) an all-sky accurate interstellar model and a spectrum of the H 1 7 -ray emissivity 


were presented, followed by various publications related to the interpretation of an apparent excess 
observed around 1 GeV in the EGRET data. 

The Large Area Telescope (LAT) is the main instrument onboard the Fermi satellite. It collects 
7 rays in the energy range 20 MeV to greater than 300 GeV. The LAT has a wide field of view 
of about 2.4 sr, an on-axis effective area of 8000 cm 2 and a 68 % containment of the point-spread 
function (PSF) of 0.8° at 1 GeV. For most of the first four years of the mission it was operated in 


all-sky survey mode and imaged the sky every two orbits in 3 hr. Abdo et al. (2009b) used LAT 


photons in a mid-latitude region in the third quadrant lacking known molecular clouds to extract 
the 7 -ray emissivity per H I in the neighborhood of the solar system for energies from 100 MeV 
to 9 GeV. They masked the bright 7 -ray point sources, and correlated for each energy bin the 
H I column density (-/V(H 1 ) ) with the LAT counts after subtracting the contribution from IC. 
With the post-launch response function Pass 6_V3 DIFFUSE, they measured an integrated 7 -ray 
emissivity of (1.63±0.05)xl0 -26 ph s - 1 sr - 1 atom -1 above 100 MeV and found that the differential 
emissivity agrees with calculations based on CR spectra consistent with those directly measured. 
The EGRET apparent excess not confirmed by Fermi (Ab do et al.| 2009a) was of likely instrumental 
origin ( Stecker et al.| 2008). Delahaye et al. (2011) compared those emissivities with predictions 
from various proton and Helium local interstellar spectra (LIS) and a 7 -ray production cross-section 


of Huang et al. (2007). They found good agreement with the LAT measurements but neglected the 


contribution from electrons. Emissivity was also derived from Fermi -LAT observations in specific 


regions like Cassiopeia and Cepheus (Abdo et al. 2010), the third Galactic quadrant (Ackermann 















































3 


et al. 2011), Orion (Ackermann et al. 2012e), the Cygnus region (Ackermann et al. 2012a) and 


Chamaeleon, R Coronae Australis and Cepheus and Polaris flare regions (Ackermann et al. 2012c 


2013). 


The objective of the present work is to derive the spectra of local interstellar CRs from Fermi- 
LAT observations independent of solar modulation which affects direct measurements in the he¬ 


liosphere. We have extended the work of Abdo et al. (2009b) to the whole longitude range and 


to latitude 10° < | 6 | < 70° to derive from 4 years of LAT data the local H I emissivity using a 
template fitting method. This method is based on the fitting of the Fermi-LAT counts map by a 
linear combination of templates spatially correlated with predicted production sites of 7 rays. We 
compare those emissivities with the ones computed from heliospheric CR observations and 7 -ray 
production cross-sections. We derive an electron-positron, proton and Helium LIS compatible with 
the H 1 emissivity using the proton and Helium fluxes detected with PAMELA (Adriani et al. 2011) 


as well as the electron plus positron fluxes measured by Fermi-LAT (Ackermann et al. 2010). 


ISM Neutral Gas Census 


Most 7 rays in the LAT energy range produced in the ISM come from CR protons imping¬ 
ing on atomic and molecular hydrogen and, to a lesser extent for the energies considered here, 
bremsstrahlung radiation when electrons and positrons are deflected by the gas nuclei. The num¬ 
ber of 7 rays produced is then expected to be proportional to the neutral hydrogen column density. 
CRs also interact with the Helium and heavier-element gases of the ISM, here we assume that 
they are uniformly mixed with the hydrogen. We used the all-sky Leiden-Argentine-Bonn (LAB) 
(Kalberla et al.| 2005) 21-cm line radiation composite survey to obtain the atomic hydrogen 7V(H 1 ) 
column densities with the assumption of a uniform spin temperature (Tg). We used the J=l—>0 
line of carbon 12 monoxide (CO) at 2.6-mm as a tracer for molecular hydrogen H 2 . The velocity 


integrated CO intensity VP(CO) was obtained from the Center for Astrophysics compilation (Dame 
et al. 2001 ) as well as new observations extending the coverage in the northern hemisphere with the 
same telescope (T. Dame 2011, private communication). In our study, we did not assume a value 
for the molecular hydrogen to CO conversion factor Xco = JV(H 2 )/W(CO); instead, we obtained 
it from the fitting procedure described in Section [ 6 j The ’’Dark Neutral Medium” (DNM) gas 
( |Grenier et al. 2005b; |Ade et al.|[2011[ |Paradis et al.|[2012[ ) , found in excess over estimations based 
on iV(H 1 ) and IT(CO), was accounted for by including in our model a template for total dust 
column-density since neutral gas and dust column-density are to the first order linearly related 
( |Boulanger et al. 1996). We used the dust optical depth map of Schlegel et al. (1998) to trace the 
total dust column density. This map is based on the IRAS/ISSA map at 100 /mi with the infrared 
intensity corrected for dust temperature based on the ratio of COBE/DIRBE at 100 and 240 fim. 
It was normalized as E(B-V) reddening expressed in mag. We modeled the dust optical depth map 
with a linear combination of 1V(H 1 ), VE(CO) and isotropic emission maps with coefficients free to 
vary. We fitted this simple model to the dust column-density map and obtained a residual dust 
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map corresponding to dust not linearly correlated with the N (H i) and IF(CO) column-densities. 
The final DNM map was obtained after applying to this residual map an absolute magnitude cut 
at 0.04 mag to remove diffuse high latitude dust residual clouds in low density regions where dust 


optical depth and hydrogen column density depart from portionality (see Figure 6 of Ade et al 


( 2011 )). 


To allow for a Galactocentric gradient of CR flux in the Galaxy, we partitioned the IV(H i) 
column-densities and W(CO) intensities into several Galactocentric annuli. The present work 
interprets the emissivities of the local H I annuli corresponding to Galactocentric distances from 
about 8 to 9.5 kpc (the IAU recommended value for the Sun Galactocentric distance is 8.5 kpc). 
Restricting the analysis to absolute latitude above 10° insured that most 7 rays originated from 
collisions between CR and hydrogen located at distances within few hundreds of parsecs from the 
sun. A more precise description of the gas annuli used in this work is given in Ackermann et al 


( 2012 d). 


3. Ionized Hydrogen 


Hydrogen also exists in form of H II regions embedded in a pervasive diffuse plasma layer 
(Reynolds 1991; Ferriere| 2001 Gaensler et al.| 2008 Haffner et al. 2009). This layer called the 
warm ionized medium (WIM), or sometimes “Reynolds Layer”, represents nearly 90% of the H + 
with a scale height of 8 times the scale height of the neutral hydrogen (Gaensler et al. 2008). First 
introduced to interpret the low-frequency absorption of Galactic synchrotron, it was later confirmed 
based on the dispersion measures of pulsars and detected through optical recombination line emis¬ 
sion (mainly Ha) and Coulomb interactions of free electrons with ions (thermal bremsstrahlung or 
free-free radiation) at microwave frequency. 

The total mass of H + in the Galaxy is estimated from the combination of dispersion measures 
of pulsars and models of the density distribution of Galactic free electrons. The model NE2001 de¬ 


veloped by Cordes & Lazio (2002) is commonly used in pulsar studies, Schnitzeler (2012) compared 


dispersion measures to predictions from eight electron density models. The conclusion was that a 
modified version of NE2001 best reproduces the observations. In NE2001 the WIM is modeled by 
a thick disk with a squared hyperbolic secant dependence on the height above the Galactic disk. 
Additionally the model includes a thin annular disk, spiral arms, local bubbles and individual re¬ 
gions of high and low electron density (clumps and voids). The model parameters are tuned to best 
reproduce 1143 dispersion measures of pulsars with known distances. Since the Helium is neutral 
in the WIM and since its contribution to the hot ionized medium electron density is limited, we cal¬ 
culated the ionized hydrogen density IV(H + ) directly as the free electron density given by NE2001. 
We obtained from NE2001 that N(H + )/N(H 1 ) equals 0.17 in the Galactic plane | 6 | < 10° and 


0.3 outside this region. Alternatively, He et al. (2013) studied 68 radio pulsars detected at X-ray 


energies and compared the free electron column density given by dispersion measures to 1V(H 1 ) 
along the line of sight traced by X-ray extinction, they obtained a ratio in the range 0.07 to 0.14. 
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As we describe below, we fitted various 7 -ray production site templates and point sources to the 
LAT data (see Sectionjfi]). We added to the fit a template of the WIM based on NE2001 predictions 
with a normalization factor free to vary, this component did not improve the model likelihood. We 
did not observe any improvement when applying latitude cuts or by removing individual H + clumps 


that seemed over predicted. Paladini et al. (2007) also observed a dust emission fit worsen with 


an H + template extracted from NE2001. This probably comes from the simplified model used in 


NE2001, Figure 7 of Sun et al. (2008) shows that NE2001 does not reproduce the structures of the 


free-free emission template from Wilkinson Microwave Anisotropy Probe (WMAP) observations. 


Based on Gaensler et al. (2008) we also created an H + template based on an exponential scale- 


height for H + leading to a cosecant IV (H + ) template map which also did not improve the model 
likelihood in the fit to the LAT observations. 

While the dispersion measure provides the column density of free electrons, the free-free and 
Ha emissions, which involve two-particle processes, are proportional to the integral of the square 
of the electron density along the line of sight (emission measure). Ha emission suffers from dust 


absorption particularly for absolute latitude below 15° (Dickinson et al. 2003) and a substantial 


fraction of the diffuse high-latitude Ha intensity is the result of scattering by interstellar dust of Ha 
originating elsewhere in the Galaxy (Witt et al. 2010| Seon & Witt 2012). The free-free emission 
does not suffer from those flaws but microwave surveys also include synchrotron radiations. We 
used the free-free intensity map at frequency v= 22.7 GHz (Iff) extracted from 9-years of WMAP 


observations with the maximum entropy method (Bennett et al. 2013) as a proxy for the 7 -ray 
emission correlated with H + . 

The conversion from emission measure to column density is done by assuming an effective 
electron density n e f f such that the integral of the square of the electron density n e along the line of 
sight is f rigds = n e ff f n e ds. Under this approximation, the free-free emission intensity in Jy sr -1 
is Iff = 2.5 X 10 6 T e ~ 0 35 i/~ 01 n e ff f n e ds where the electron temperature T e is in K, the frequency 
v in GHz and the line of sight distance in kpc (Sodroski et al. 1997). We can then express the H + 
column density as: lV(H + )=/’ n e ds = 1.2 x 10 15 T e 035 v 01 n e ff~ 1 Iff. We used !V(H + ) in Equation 
[l] adopting an electron temperature T e =8000 K (Sodroski et al. 1989). 


4. Inverse Compton and Large Scale Structures 

The inverse Compton scattering of the interstellar radiation field by electrons and positrons 


represents a non-negligible fraction of the 7 -ray interstellar emission (Ackermann et al. 2012d) 


There is no simple template for the IC emission so we relied for its morphology on the program 
GALPROFp] ( |Strong et alT 2007, Vladimirov et al.||2011 ) which predicts the CR density by solving 
numerically the transport equation of charged particles, starting from a given source distribution 


1 http://galprop.st anford.edu 
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and propagation parameters. This density is combined by GALPROP with an interstellar radiation 
field model (Porter et ah 2008) into a predicted IC intensity map Iic p - We used an IC morphology 
predicted by one of the diffusive reacceleration models described in Ackermann et al. (2012d) with 


the identification 5 y^6' R 30 r 150 <:: 2 (CR source distribution proportional to the radial distribution 


of pulsars in the Galaxy given by Yusifov & Kiigiik (2004), Galactic halo size and radial boundary 


equal to 6 kpc and 30 kpc). In Section [ 8 ] we also used two different IC templates also derived 
from GALPROP diffusive reacceleration models but with different input parameters designated 
54_zlOG4c5rS (Gaussian source distribution with a halo height of 10 kpc) and 54_77Xvarh7S (source 
distribution intermediate between the radial distribution of SNRs and that of pulsars, halo height 
of 4 kpc). During the fitting procedure, we left the Iic p normalization free in each energy bin to 
allow some freedom for possible variations of the CR distribution and Galactic interstellar radiation 
field from the GALPROP model. 

GALPROP does not include local structures like the Loop I giant radio continuum loop. 
Gamma rays likely coming from IC emission of CR electrons trapped in Loop I have been reported 
( |Grenier et al. 2005a; Casandjian Sz Grenier 2009 Ackermann et al. |2012d ). We modeled this 
emission using the 408 MHz radio-continuum emission survey ( iLoopi ) of Haslam et al. (1982). We 


subtracted the radio point sources from this map and created a dedicated template from a selected 
region around Loop I. Since the 408 MHz radio-continuum emission depends on the magnetic field 
while the 7 emission does not, we cannot rely only on the radio template for our analysis. That 
is, we cannot assume that the 7 -ray intensity is proportional to the radio continuum intensity. We 
extended the radio template of Loop 1 with a patch of uniform intensity ( I pa tche ) with corners located 
at ( l,b ) = (350°, 0°), (310°, 25°), (300°, 17°), (332°,—12°). This shape approximately reproduces the 
extension observed map of residuals between the LAT observations and the model when only the 
408 MHz radio-continuum template was used, this patch possibly compensates for the difference 
between the radio and 7 -ray emission of Loop 1 in a region where the magnetic field might be 
weaker. 


Additionally we have incorporated in our model a template for the Fermi Bubbles (Su et al. 


2010; Ackermann et al. 2014). There is no accurate a priori template for the 7 -ray emission of those 
large structures, but their morphologies are distinct from those of the gas or of the GALPROP IC 
templates and their spectra are harder. We derived their shapes from an iterative procedure for 
energies integrated above 4 GeV, starting from a fit using a uniform intensity template with the 
same shape as that of the Bubbles. We added to this template the fit residuals and verified that 
the template derived above 4 GeV {IBubbles) gives a flat residual map at lower energies. 


5. Data Selection 


Our analysis is based on Pass 7 reprocessed data (P7REP) Clean class events for the first 
4 years of the mission. We excluded from this time interval periods during intense transients 


(fluent GRBs, bright AGN or solar flares) occurred (Ballet & Burnett 2013). We used the Clean 
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selection that has a reduced residual charged-particle background compared to the Source selection 
(Ackermann et al. 2012b). The event selection included 55 million LAT counts which were binned 
into 26 energy intervals logarithmically spaced from 50 MeV to 50 GeV. For energies below 50 MeV, 
at which the Pass 7 effective area is approximately 20 times lower than at 1 GeV, we did not obtain 
a good fit stability. The lower statistics, a higher Earth limb contamination and a broadening of the 
LAT point-spread function (PSF) led to confusion between point sources and Galactic interstellar 
emission. The highest four energy bins were combined into two bins. In each energy interval 
the high statistics allowed a good separation of photons correlated with H I from those with other 
origins for absolute latitude above 10°. Figure[l](a) shows the total Fermi -LAT 4-year accumulated 
number of counts above 170 MeV ( N obs ) for the latitude range considered here. 


The residual intensity from Earth limb photons entering the LAT at large zenith angles ( Iu m b ) 
cannot be neglected in the lowest energy bins. It is structured around the celestial poles and can 
bias the fit. We derived a template for this residual contamination of limb photons by subtracting 
a counts map derived with a zenith angle cut at 100° to one restricted at angles above 80°. We 
added this contribution to the model as a template free to vary. 


6. Fitting Procedure 


The template method, based on fitting the 7 -ray observations with a linear combination of 
maps spatially correlated with production sites of 7 rays (templates), has been widely used in 
diffuse 7 -ray astronomy (e.g., Lebrun & Paul 1979 Strong et al. 1988 |Digel et al. 1999 


Grenier 


et al. 2005b: Abdo et al. 2010). The fitting procedure is similar to the one used to build the Fermi- 


LAT interstellar emission model 0 We modeled the predicted spatial distribution of counts {N pre d) 
detected by Fermi-LAT as: 


N pre d(E,l,b) = qi(E)Ii(l,b) + Nic{E)Iic p {E,l,b) + Ni S0 (E)Ii 30 + Nl oopl (-^) 1 LoopI (J< ? fr) 

i=gas 

+ N pa tch(E)I p atch(l, b) + NBubbles(E)lBubbles(l , b) + N pti ( E)8(l , b, l) 

i=point src 

+ ^ ] Nexti^E^Iid^b) + I S un_moon(E,l,b) + Nlimb(E)Ilimb(l, b) 

i=extend src 


(i) 


where E is the energy and l,b are Galactic sky coordinates. In Equation [l] the term “gas” refers 
to IV(H 1 ), VF(CO), IV(H + ) and the DNM column densities. They are multiplied by the respective 
hydrogen 7 -ray emissivities q % . The factors N IC , N iso , N LoopIl N patchl N Bumes and N Umb represent 
the renormalization factors associated with the predicted IC, the isotropic emission represented 


“http://fermi.gsfc.nasa.gov/ssc/data/access/lat/Model_details/FSSC_model_diffus_reprocessed_vl2. 
pdf 
























Fig. 1.— a: Fermi -LAT counts per pixel integrated above 170 MeV together with predictions cal¬ 
culated from Equation [l] for 7 rays originating from local H I (b) and from remaining contributions 
(c). When the contributions from (b) and (c) are subtracted from the total counts map (a) we 
obtain the residual plot (d) displayed in units of standard deviations. We used the same scale for 
maps (a-c). We masked the absolute latitudes below 10° and above 70° corresponding to region 
outside our region of interest. In this figure the maps are displayed in Mollweide projection in 
Galactic coordinates, we applied a smoothing with a gaussian symmetric beam of 2° (FWHM). 
The pixel size for all the maps is 0.25°. 
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by a uniform intensity template, the Loop I radio emission and patch, the bubbles and the limb 
template maps. The normalization of the template^] corresponding to the extended emission of 
the two Magellanic Clouds and of Centaurus A are written N ex t z . The 2179 point sources listed 
in a preliminary version of the 3FGL catalog (Ballet & Burnett 2013) are also fitted individually 
in Equation [lj it is represented by a number of counts N pbi and a Dirac 5 function at the source 
position. In Equation [l] we denote by I the counts map resulting from the convolution with the 
LAT PSF of the product of an intensity map I by the instrument exposure and pixel solid angle. 
We obtained the gas emissivities and normalization factors from a binned maximum-likelihood fit 


with Poisson statistics to the Fermi -LAT counts map in the 26 energy intervals (Mattox et al. 


1996). We used the code MINUIT^to calculate the likelihood maxima as well as the parameters 
error. The intensity from the Sun and the Moon (Abdo et al. 2011, 2012), I S un_moom is kept fixed 
in the analysis. We restricted the fit to latitudes 10° < |6| < 70° but included in Equation [I] lower- 
latitude sources and gas templates to account for the spread of the bright Galactic ridge above 10° 
at low energy. We excluded regions with absolute latitudes above 70° where the isotropic emission 
dominates. We binned all the maps in HEALPrx^] with a Nside of 256 which corresponds to a bin 
size of about 0?25 x 0?25. In each energy interval N IC , N iso , N LoopI , N patch , N BubUes , N pti , 
Nexti, and N[ tmb were left free to vary. 


7. Results 

We fitted the model given by Equation [T] to the Fermi-LAT counts in each of the 26 energy 
bins. We derived for each template and source an intensity that we averaged over the region studied 
here and plotted in Figure [2] with vertical error bars obtained from the likelihood maximization. 
In this graph we also show the intensity calculated from the LAT counts for the same region with 
the corresponding statistical uncertainties. Horizontal error bars represent the bin energy widths. 
To simplify the graph we combined the various gas contributions into a single spectrum called 
“ 7 r°, bremsstrahlung”; for this spectrum we combined quadratically the uncertainties of individual 
contributions. We observe a very good agreement between the LAT observations and the sum of 
the template model intensities. We note that 7 rays originating from the interaction between the 
CRs and the gas dominate for energies between 200 MeV and 10 GeV. 

We study and interpret in this paper the contribution in Equation [l] from the atomic hydrogen 
template: qHiN(R 1 ). In Figure [l]^ we show this contribution, in counts per pixel integrated 
above 170 MeV. In Figure [lj: we show the complementary contribution from all the other terms of 
Equation [l N W ed ~ QhiN(H. 1 ). When subtracting all the contributions to the LAT counts map we 


3 http://fermi. gsfc.nasa.gov/ssc/data/access/lat/2yr .catalog/ 

J h ;tps://wwwasdoc.web.cern.ch/wwwasdoc/minuit/minmain.html 
‘http://healpix.sourceforge.net 
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Fig. 2.— Intensity spectrum for the LAT counts and individual contributions of Equation [I] for 
latitudes 10° < |6| < 70°. The uncertainties associated with the LAT intensity are only statisti¬ 
cal. The vertical error bars associated with each contribution represent the uncertainty obtained 
from the fit. For the most part they are smaller than the symbol size. The contributions asso¬ 
ciated with hydrogen in atomic, molecular, ionized, or DNM phase are merged into a single a 7r°, 
bremsstrahlung” contribution. 
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obtained a flat residual map (N 0 b s — Np re d)/ y/N pre d (Figure [ljl) which indicates that the main LAT 
7-ray sources are included in Equation [l] and then that the local H I template does not compensate 
for any un-modelled emission. 

In each energy bin we calculated the map of residual LAT counts without including the local 
H i template in the model: N 0 b SHI = N a b s — ( N pre d — qmN (H i)), and divided by the exposure 
map and pixel solid-angle to obtain the intensity I 0 bs HI ■ This intensity represents the intensity of 7 
rays measured by the LAT not associated by the fit to any model components other than the local 
H 1. In Figure [3] we show the pixel-by-pixel correlation between I 0 bs HI and the predicted intensity 
QhiN {H i). In this figure we have divided our total energy range in four intervals, and averaged 
the intensity within 256 neighboring pixels (1024 pixels for the highest energy range) to lower the 
Poisson fluctuation. We observe a very good linear correlation between 7 rays and atomic hydrogen 
column density. 


The instrument energy dispersion introduces a bias between the number of reconstructed counts 
in a given energy range and the number of incoming photons of corresponding energies. Above 
200 MeV, for LAT event selection P7REP Clean, the combination of the incoming photon spectrum 
and the energy dependence of the LAT acceptance reduces the effect of energy dispersion; counts 
that are dispersed by an amount larger than half an energy bin width are equally distributed in 
bins with energies higher and lower than that of their corresponding photons. Below 200 MeV the 
sharp fall-off of the 7r° decay spectrum combined with the decrease of the effective area of the LAT 
breaks this equilibrium and results in an increase of counts in the lowest energy bins. We used the 
LAT Science Tooj^] gtrspgen to estimate this effect for photon energies from 50 MeV to 50 GeV and 
observed that a correction term for the bias introduced by the energy dispersion follows the simple 
formula l + 7500il ,_2 ' 43 . It leads to an increase of counts of 38% in our lowest energy range spanning 
from 50 to 64 MeV and of 21% for the one from 64 to 82 MeV. The effect of the energy dispersion 
is less than 2% for the range spanning from 172 to 220 MeV. We proceeded iteratively to account 
for the dispersion, we first used a predicted diffuse spectrum to calculate an approximate energy 
dispersion correction. We then applied this correction to derive the H 1 emissivity, and used the 
measured spectrum for a more accurate dispersion correction. We verified that no extra iteration 
was necessary. 


Taking into account this correction for the energy resolution we plotted in Figure [4] the local 
7-ray differential emissivity per H 1, obtained by dividing the emissivities qm by the energy band 
width. We represented the statistical uncertainties (see Section 8) with vertical error bars and 
instrumental systematic uncertainties with brackets. We positioned the points at energies corre¬ 
sponding to the logarithmic mid-points of the bands and scaled the emissivities by the square of 
those energies. We include the differential emissivities from previous work from COS-B, EGRET, 
and previous Fermi -LAT results. The agreement is good except for the EGRET emissivities above 
1 GeV as studied by Abdo et al. (2009a). The differential emissivities multiplied by the square of 


6 Available from http://fermi.gsfc.nasa.gov/ssc/ 
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Fig. 3.— Pixel by pixel correlation in 4 energy bands between the 7-ray intensity and H 1 column 
density. We multiplied the LAB survey column density by the emissivity per H 1 and plotted this 
value against the LAT counts not associated by the fit to any model components other than the 
local H I divided by the exposure to yield corresponding intensities. We selected pixels located 
within an absolute latitude range of 10° to 70°. Each point corresponds to a pixel size of 3.7° (7.3° 
for the highest energy range). 
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the band mid-point energy are provided in Table [lj 

We investigated possible variation of the emissivity spectrum across the sky. We divided our 
region of interest in four quadrants and restricted the fit of Equation[l]to the quadrant North-West: 
10° < b < 70° and 180° < l < 360°, South-West: -70° < b < -10° and 180° < l < 360° North- 
East: 10° < b < 70° and 0° < l < 180°, and South-East: —70° < b < —10° and 0° < l < 180°. We 
show in Figure[5]the emissivities in the four quadrants divided by the emissivity obtained previously 
in the whole region. We obtain a good agreement between the emissivities in the quadrants, 
indicating no significant local CR asymmetry across the sky. 

We verified that the 7-ray emissivities measured in the local atomic hydrogen are compatible 
with the ones extracted from the molecular hydrogen, DNM gas and ionized hydrogen. Similarly 
to what we did for H 1 we obtained from the fit of Equation [I] the emissivities associated to W(CO) 
( qco ) and to the DNM template ( qDNM ) and derived the molecular hydrogen-to-CO conversion 
factor Xco = Qco/^Qhi and the dust-to-gas ratio Xoust = Qdnm/qhi- Restricted to the energy 
range above 200 MeV, where the correlation between the CO template and point sources is weaker, 
we obtained Xco = (0.902 ± 0.007) x 10 20 cm -2 (K km s -1 ) -1 and Xoust = (41.4 ± 0.3) x 10 20 
cm -2 mag -1 . The values quoted here correspond to clouds and dust within the absolute latitude 
range 10° to 70°. The uncertainties do not account for systematics (see Section [8]). Figure [6] 
shows the correlation between the 7-ray emissivity in the atomic and the molecular phase (a) and 
DNM (b) obtained from qco and qoNM scaled with the ratios Xco and Xoust ■ As in Figure 
[4] we scaled the emissivities by the square of the energy. We observe a good linear correlation 
(Figure [6^,-b) indicating that the spectral shape of the hydrogen emissivity is similar when the 
hydrogen is atomic, a molecule and in the gas that is traced solely by dust. To verify the method, 
we proceeded similarly with the isotropic template by scaling lVj so by a factor Xi so to match the 
range of the emissivity per H I. As expected there is no linear correlation between the isotropic 
normalization and the atomic hydrogen emissivity (Figure [Gji) . 

In Figure [6^ we show for the first time the linear correlation between 7-ray emissivities mea¬ 
sured in H 1 and in H + ( qh +)- As explained in Section[3]we converted the ionized hydrogen free-free 
emission measure to H + column density. This conversion depends on the the electron effective den¬ 
sity parameter n e ff that we adjusted to match the emissivity in H + to the one in H I. We obtained 
the best fit for n e //=4.6 cm -3 which is in the expected range of expected (Sodroski et al. 1989). 
That corresponds to !V(H + )= 1.3 x 10 20 7// where Iff is expressed in mK and N(H + ) in atoms 
cm -2 and to an average ratio N(H + )/N(K 1) equal to 0.015. This proportion of H + to H I atoms is 
notably less than that estimated from dispersion measures (see Section [3]). The fit may be biased 
by the strongly emitting H 11 regions that are incorrectly reconstructed by our average electron 
temperature. The remaining 7-ray emission correlated with H + might then be associated to other 
smoothly varying templates of Equation [T] like the IC or the isotropic distribution. It is also possible 
that the H + is traced together with the DNM by the dust or that it follows the H 1 distribution. In 
this case the 7-ray emissivity per H 1 measured in this work would be overestimated by the amount 
of H + mixed in H 1. The large uncertainties at low energies of the hydrogen emissivities measured 


Sodroski et al. 1989 
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Fig. 4.— Differential local 7-ray emissivity per H I versus 7-ray energy. The vertical error bars 
correspond to statistical uncertainties while systematics are represented by square brackets, the 
horizontal bars correspond to the bins energy width. The differential emissivity is compared to 
previous works from COS-B (Strong et al. 1988), EGRET ( |Strong fe Mattox 1996), and to early 
Ferrra-LAT results (lAbdo et al. 2009b). 






















Table 1. Local H I differential emissivity 


Energy band 
MeV 

E 2 (qn 

r l zb stat zb sys ) 

( a ) 

E 2 ( qco ± stat ) 

(b) 

E 2 {qDNM ± stat ) 

(c) 

E 2 ( q H + =b stat ) 

(d) 

NiC p 

zb stat 

E 2 ( N iso ± stat ) 

(e) 

50-64 

0.47 

± 

0.02 

± 

0.09 

1.27 

± 0.10 

0.00 

± 0.01 

0.0 

± 0.6 

0.60 

± 0.01 

2.03 

± 0.02 

64-82 

0.66 

± 

0.02 

± 

0.11 

1.23 

± 0.09 

0.06 

± 0.01 

0.0 

± 0.5 

0.71 

± 0.01 

1.85 

± 0.01 

82-105 

0.82 

± 

0.02 

± 

0.11 

1.50 

± 0.09 

0.25 

± 0.02 

0.0 

± 0.4 

0.84 

± 0.01 

1.73 

± 0.01 

105-134 

1.09 

± 

0.01 

± 

0.13 

1.85 

± 0.08 

0.36 

± 0.02 

0.8 

± 0.3 

0.90 

± 0.01 

1.61 

± 0.01 

134-172 

1.36 

± 

0.01 

± 

0.14 

2.45 

± 0.08 

0.49 

± 0.02 

1.0 

± 0.4 

1.04 

± 0.02 

1.59 

± 0.01 

172-220 

1.71 

± 

0.01 

± 

0.15 

2.84 

± 0.07 

0.66 

± 0.02 

1.9 

± 0.4 

1.18 

± 0.02 

1.54 

± 0.01 

220-281 

1.99 

± 

0.01 

± 

0.16 

3.44 

± 0.07 

0.75 

± 0.02 

2.9 

± 0.3 

1.29 

± 0.02 

1.44 

± 0.01 

281-360 

2.22 

± 

0.01 

± 

0.16 

3.98 

± 0.07 

0.90 

± 0.02 

2.9 

± 0.3 

1.43 

± 0.02 

1.32 

± 0.01 

360-461 

2.37 

± 

0.01 

± 

0.15 

4.28 

± 0.07 

0.96 

± 0.02 

2.4 

± 0.3 

1.60 

± 0.02 

1.14 

± 0.01 

461-589 

2.45 

± 

0.01 

± 

0.13 

4.32 

± 0.07 

1.03 

± 0.02 

3.6 

± 0.3 

1.72 

± 0.02 

1.01 

± 0.01 

589-754 

2.46 

± 

0.01 

± 

0.13 

4.39 

± 0.07 

1.03 

± 0.02 

3.1 

± 0.3 

1.79 

± 0.03 

0.90 

± 0.01 

754-965 

2.39 

± 

0.01 

± 

0.14 

4.33 

± 0.07 

1.05 

± 0.02 

3.7 

± 0.3 

1.90 

± 0.03 

0.75 

± 0.01 

965-1235 

2.29 

± 

0.02 

± 

0.14 

4.20 

± 0.08 

0.96 

± 0.02 

2.8 

± 0.3 

1.95 

± 0.04 

0.64 

± 0.01 

1235-1581 

2.12 

± 

0.02 

± 

0.14 

3.98 

± 0.08 

0.87 

± 0.02 

2.9 

± 0.4 

2.10 

± 0.04 

0.55 

± 0.01 

1581-2024 

1.97 

± 

0.02 

± 

0.14 

3.65 

± 0.08 

0.82 

± 0.02 

2.1 

± 0.4 

2.14 

± 0.05 

0.52 

± 0.01 

2024-2590 

1.70 

± 

0.02 

± 

0.13 

3.03 

± 0.08 

0.73 

± 0.02 

2.2 

± 0.4 

2.18 

± 0.06 

0.45 

± 0.01 

2590-3314 

1.42 

± 

0.02 

± 

0.13 

2.67 

± 0.08 

0.61 

± 0.02 

2.0 

± 0.4 

2.12 

± 0.05 

0.41 

± 0.01 

3314-4242 

1.20 

± 

0.02 

± 

0.12 

2.22 

± 0.08 

0.54 

± 0.02 

1.2 

± 0.4 

2.26 

± 0.06 

0.36 

± 0.01 

4242-5429 

1.02 

± 

0.02 

± 

0.11 

1.85 

± 0.09 

0.43 

± 0.02 

1.6 

± 0.4 

2.25 

± 0.06 

0.34 

± 0.01 

5429-6947 

0.83 

± 

0.02 

± 

0.10 

1.64 

± 0.09 

0.35 

± 0.02 

0.5 

± 0.4 

2.16 

± 0.07 

0.36 

± 0.01 

6947-8891 

0.67 

± 

0.02 

± 

0.09 

1.26 

± 0.09 

0.31 

± 0.02 

1.8 

± 0.5 

2.06 

± 0.08 

0.34 

± 0.01 

8891-11379 

0.55 

± 

0.02 

± 

0.07 

1.03 

± 0.09 

0.22 

± 0.02 

0.1 

± 0.4 

2.04 

± 0.09 

0.32 

± 0.01 

11379-14563 

0.45 

± 

0.02 

± 

0.06 

0.87 

± 0.10 

0.17 

± 0.02 

0.7 

± 0.5 

1.68 

± 0.15 

0.34 

± 0.01 

14563-18638 

0.36 

± 

0.02 

± 

0.05 

0.86 

± 0.10 

0.12 

± 0.03 

0.3 

± 0.5 

2.24 

± 0.15 

0.26 

± 0.01 

18638-30527 

0.31 

± 

0.02 

± 

0.03 

0.56 

± 0.08 

0.12 

± 0.02 

0.8 

± 0.4 

1.83 

± 0.17 

0.23 

± 0.01 

30527-50000 

0.20 

± 

0.02 

± 

0.02 

0.45 

± 0.09 

0.07 

± 0.02 

0.6 

± 0.5 

2.28 

± 0.16 

0.16 

± 0.01 

Note . — ( a ): 10 _24 MeV 2 s " 1 
MeV ” 1 mag " 1 , ( d ): 10" 24 MeV 2 

sr - 1 MeV - 1 , ( b ): 10- 4 MeV 2 s " 1 cm " 2 sr " 1 MeV " 1 (K km S - : 
s _1 sr -1 MeV -1 , ( e ): 10 -3 MeV 2 s _1 cm -2 sr -1 MeV -1 

L " 1 , ( c ): lQ — 2 MeV 2 s " 1 

— 9 — 1 

cm sr 
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Fig. 5.— Ratio of the emissivity in four sky quadrants to the emissivity in the whole region analyzed 
versus 7-ray energy. Only statistical uncertainties are shown. See the text for the definitions of the 
quadrants. 
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in H + prevented the study of screening effect in the Bremsstrahlung production function (Hunter 


et al. 1997). 


8. Systematic Errors 


We investigated three major sources of systematic errors: the H I spin temperature (Tg), the 
modeling of the IC, and the absolute determination of the LAT effective area. 

When deriving IV(H i) from the LAB radio survey we corrected for the 21-crn H I line opacity 
which depends on its Tg. Since CRs interact with H I independently of Tg and since the Galaxy 
is optically thin to 7 rays at energies considered here we investigated the uncertainty related to Tg 
by comparing the LAT counts map with models based on 1V(H 1) derived with various Tg values. 
The thinness of the local H 1 does not make JV(H 1) sensitive enough to variation of Tg therefore 
we selected a region toward the Galactic anticenter (90° < l < 270°, \b\ < 70°) to obtain the H 1 Tg 
that gives the best fit to the data. We investigated seven values of Tg from 90 K to 400 K. For each 
temperature, we derived an H 1 column density and a DNM map. We used those maps in a fit to the 
Fermi -LAT counts in the anticenter region. We observed that the likelihood was maximum for the 
model derived with Ts=140 K. We adopted this temperature to derive N(H 1) for the whole Galaxy. 
We used the values 125 K and 170 K to estimate the uncertainty in the emissivity related to the 
Tg. This uncertainty range is larger than the one obtained by varying the maximum likelihood by 
1 i.e., the statistical uncertainty in the fit parameter. This range corresponds to values of Tg above 
which a noticeable discrepancy is observed in the fit residual map. We used in Equation [l]lV(H 1) 
derived with a Tg of 125, 140 and 170 K and compared in each energy bin the resulting emissivities. 
We obtained a systematic uncertainty related to Ts of 2% corresponding to the average variation 
of those emissivities. Studies of H 1 absorption against background radio sources have shown that 


Tg is not uniform in the multi-phase ISM (Heiles &; TYoland 2003 Kanekar et al. 2011), but we 


assumed here a uniform Tg and corrected this approximation by representing in the model maps 
from positive and negative residuals of the total dust column density. 


As shown by Ackermann et al. (2012d) there is no unique set of GALPROP propagation 


parameters that can model the Galactic diffuse emission in the Fermi-LAT sky. Since we rely on 
GALPROP for the morphology of the IC map, we investigated the effect of three sets of GALPROP 
parameters labelled 54_zlOG4c5rS, 54_77Xvarh7S and S Y Z 6 R 30 T 150 C 2 (Ackermann et al.(2012d). 
This method does not account for all the uncertainties related to the Galprop IC map but since 
the spatial structures of IV(H 1) allow for a good separation from the smooth morphology of the 
IC map, we do not need to input into our fit an accurate IC template. As an illustration, when 
we do not include into the fit any IC template, the HI emissivity varies by only 20%. The effect 
of using the three different IC models is limited, from 3% at low energy to 1% at higher energy. 
We used recent IC distribution predictions obtained with the version S Y Z 6 R 30 T 150 C 2 for the final 
7-ray fit and added in quadrature to the systematic uncertainties the variation obtained from the 
other parameter sets. 
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Fig. 6.— Comparison between qni and (a) qco , (b) qDNM , (c) (/#+, and (d) the isotropic nor¬ 
malization. The emissivities are scaled by the square of the energy at which the emissivity was 
measured. The open circle symbol represents the emissivities derived for energies below 200 MeV. 
The emissivities measured in the molecular, ionized and DNM are scaled to the emissivity per 
H I so that identical emissivity spectra in those different phases would be aligned on the dashed 
median line. We observe a similar spectral shape for the hydrogen emissivity in the four compo¬ 
nents: atomic, molecular, ionized, and in the DNM. As expected we do not observe any correlation 
between the isotropic normalization factor and the emissivity since no correlation is expected be¬ 
tween extragalactic 7 rays and those issued from the atomic hydrogen. Error bars correspond to 
statistical errors only. 
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Finally we used the following values for the systematic uncertainties related to LAT effective 
area: 20% at 50 MeV, 10% at 100 MeV, 5% at 560 MeV, 10% at 10 GeV and we extrapolated 
linearly those values using the logarithm of the energy. Those values are more conservative than 
the recommendecQ ones. The LAT effective area uncertainties dominate the uncertainties related 
to Ts and IC in the whole energy range. 

There are as many 7 rays originating from interstellar emission within 6° of the Galactic plane 
as of the rest of the Galaxy. At low energy, where the PSF is broad, counts corresponding to 7 
rays originating from the plane spread out to higher latitude. We incorporated into Equation [l] 
the inner gas annuli that model the Galactic plane to account for this spreading. We tested the 
stability of the fit with four cuts corresponding to absolute latitudes higher than 5°, 7°, 10° and 
15° and found no statistically significant variation. We also excluded all 7 rays detected at absolute 
latitudes higher than 70°, where the number of counts coming from collisions with H 1 is four times 
less than those of extragalactic and instrumental origins. For the final emissivity fit we applied a 
cut at 10° which allows a good template separation by the fit while excluding most of the Galactic 
plane. 

Incorrect modeling of large scale structures like Loop 1 or the Fermi Bubbles would bias the 
fit. To investigate the influence of those large structures, we restricted the fit to a region away from 
them at longitudes between 50° and 280°. We observed no noticeable difference for the local H 1 
emissivity. This test, in a region where the intensity of 7 rays from IC is much lower than from 7r° 
decay, also indicates that correlation between Ijc p and IV(H 1) maps have negligible effects on the 
measured emissivities. 


9. Interpretation 
9.1. Production cross-section 


Most 7 rays with energies between 100 MeV and 50 GeV originate from the decay of 7r° 
produced in hadronic collisions when CR protons with energies from 0.5 to 10 3 GeV interact with 
ISM nuclei. At those energies, the vast majority of 7r° are produced with a low transverse momentum 
transfer through soft interactions where the large strong force coupling constant prevents the use 
of a perturbative quantum chromodynamics approach. This regime is poorly understood and the 
7T° production cross-section is obtained mainly through phenomenological models. In this work we 
have tested three 7-ray production cross-sections: 


In Kamae et al. (2005) and Kamae et al. (2006) the authors have used the accumulated 
knowledge of high-energy accelerator physics to interpret the 7r° production cross-section for as- 
trophysical purpose. They used the event generator Pythia (Sjostrand et al. 2006) to evaluate the 


1 http://fermi.gsfc.nasa.gov/ssc/data/analysis/LAT_caveats.html 
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non-diffractive component above a p-p center of mass energy of 10 GeV where phenomenological 


fits based on the Regge theory (Donnachie & Landshoff 1992) are believed to be valid. At lower 


energy they used the scaling formula of Stephens & Badhwar (1981) with parameters based on 


Blattnig et al. (2000) but readjusted to account for the additional model components. Pythia of¬ 


fers an implementation of phenomenological QCD that effectively takes care of higher order effects 
(e.g., violation of Koba-Nielsen-Olesen scaling) that covers most of the deviation from a simple 
scaling law in the non-diffractive cross-section. Single and double diffraction were considered as 
well as resonances at 1232 MeV/c 2 and around 1600 MeV/c 2 . This model provides a good fit of 
the experimental 7 r° inclusive cross-section. 


A similar approach was used by Huang et al. (2007) where the cross-section for proton energies 


above 20 GeV was calculated using the Monte Carlo particle collisions code DPMJET-III (Roesler 


et al. 2001) that describes the nuclear interactions with the two-component Dual Parton Model 


again based for soft processes on the Regge phenomenology. At lower proton energies Huang et al. 


(2007) used the same parametric model as Kamae et al. (2006) and incorporated to their model 
resonances at 1232 MeV c -2 and 1600 MeV c -2 . The contribution of g mesons (classified with the 


direct products in Huang et al. (2007)) appears to be overestimated by DPMJET-III; we reduced 


it by a factor 8 to match the experimental cross-section measured in Aguilar-Benitez (1991). 


The approach chosen by Shibata et al. (2014) is different and based only on 7 -ray production 


cross-sections obtained in p-p collisions spanning the range of energies from the threshold pion 
production to energies of the Large Hadron Collider experiments. The cross-section is parametrized 
with an empirical formula depending on the 7 -ray energy and transverse momentum as well as four 
physical parameters that are adjusted so that the empirical cross-section and deduced quantities 


fit the experimental data (Shibata et al. 2013. Sato et al. 2012 ; Shibata et al. 2007 Suzuki et al. 


2005). 


Fermi detects 7 rays resulting not only from p-p interactions but also from the interactions of 
CR nuclei with the ISM nuclei. We used Mori (2009) to scale a pp to nucleus-nucleus cross-section. 


Table [ 2 ] gives the scaling factors of Mori (2009) obtained by comparing the 7 -ray yield from p-p, p- 


nucleus, nucleus-p and nucleus-nucleus collisions calculated with DPMJET-III. We used the values 
calculated at 10 GeV/nucleon and neglected their dependence on energy which is of the order of 
a few percent (Kachelriess et al. 2014). Table [3] lists the scaling factor for collisions induced by 


proton and Helium CRs after including an ISM relative abundance of p : He : CNO : Mg Si : Fe = 
1 : 0.096 : 1.38 X 10“ 3 : 2.11 x 10“ 4 : 3.25 x 10“ 5 (Meyer 1985; Mori 2009). From this table we 
see that: ct p jsm = 1.387 a pp and OHe.iSM = 5.12a pp . CRs heavier than Helium colliding with the 
ISM also contribute to the total 7 -ray yield. Their relative abundances were estimated at 10 GeV 
from Table 1 of Honda et al. (2004) as p : CNO : MgSi : Fe = 1 : 0.0032 : 0.0016 : 0.00037 at 
10 GeV. Using these abundances together with cross-section scaling factors of Table [2j we obtained 
a correction equivalent to 0.116<7 PP that we added to ct p jsm so that cr p +heavy,iSM = 1.503 a pp . For 
information, using a CR ratio of 0.055 between Helium and proton CR fluxes at 10 GeV, we obtain 
a total enhancement of 1.78. Here we do not use this value directly since we treated protons and 
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Helium CRs separately. 


The 7 rays detected by Fermi -LAT at energies relevant for this work are also produced when 
CR electrons and positrons of about 0.1 to 10 GeV are bremsstrahlung scattered by ISM nuclei. 
We used Gould (1969) to calculate this bremsstrahlung radiation. For this we used tabulated 01 
and 02 , and we assumed crbremssiH^) = 2 crbremss{H !)■ 


9.2. Comparison with Proton and Helium Heliospheric Spectra 


We first compared the 7 -ray emissivity per hydrogen atom measured by Fermi -LAT with the 
one calculated from CR fluxes measured directly (proton and Helium) and indirectly (electron). 
To calculate the emissivities we folded the 7 -ray production cross-sections with the proton, Helium 
and electron flux. We used the spectra of proton and Helium gathered by PAMELA between 2006 
and 2008 (Adriani et al. 2011) and ATIC-2 (Panov et al. 2006). Since PAMELA and AMS-01 
(Aguilar et al. 2002) measurements are in agreement below 10 GeV, we extended the PAMELA 
proton spectrum to energies lower than 0.4 GeV using the two lowest-energy points of AMS-01. 


The electrons that produce 7 rays at Fermi-LAT energies also spiral around the Galactic mag¬ 
netic field lines and produce radio synchrotron with frequencies from a few MHz to approximately 
ten GHz. This range corresponds to the one Webber & Higbie (2008) used to derive an electron 


LIS. We folded the bremsstrahlung cross-section with the CR electron local spectra of Webber & 


Higbie (2008) globally scaled by a factor 0.64 to match the Ferrm-LAT electron-positron spectrum 


of Ackermann et al. (2010) above 10 GeV. 


In Figure [7] we show the 7 -ray emissivity calculated with the Kamae et al. cross-section com¬ 
pared the results derived here from Fermi-LAT observations. We decomposed the total emissivity 
into contributions from bremsstrahlung and hadronic collisions. For this graph we did not solar- 
demodulate the proton or Helium spectra, but instead used directly the fluxes from PAMELA and 
ATIC-2. We used a power-law interpolation between the measurements in calculating the fold¬ 
ing. We observe a disagreement between the predicted and the measured emissivities that can be 


Table 2. 


Multiplication factors at 10 GeV/nucleon from Mori (2009) by which the p-p 


cross-section is scaled to obtain a nucleus-nucleus cross section. 


Projectile 



Target 


H 

He 

CNO 

Mg-Si 

Fe 

H 

1 . 

3.81 

11.6 

20.1 

38.9 

He 

3.68 

14.2 

42.3 

73.4 

143.1 

CNO 

11.7 

42.5 

120.6 

204.2 

386.3 

Mg-Si 

20.3 

73.2 

204.2 

343.1 

628.8 

Fe 

38.8 

142. 

384.4 

634. 

1067. 
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attributed to the solar modulation of proton and Helium CRs. 


In Figure [ 8 ] we compare our results with the emissivities calculated using the production cross- 
sections of Karnae et al, Huang et al., and Shibata et al.. We used the force-field approximation 
(Gleeson & Axford 1968) to demodulate the proton and Helium spectra with an arbitrary modu¬ 
lation potential of <f> = 500 MV. The agreement between the observed and modeled emissivity is 
greatly improved by taking into account the modulation due to the Sun. We observe a spread with 
a standard deviation of the order of 10 % between the calculated emissivities due to the differences 
in the production cross-sections. 


9.3. Derivation of the Proton, Helium, and Electron Local Interstellar Spectra 


We saw above that solar-demodulated heliospheric fluxes folded with the 7 r° cross-section 
result in 7 -ray emissivities similar to the ones we measured. Going one step further, we deduced 
the LISs for proton, Helium and electron-positrons by fitting variable LIS spectral shapes to direct 
measurements of modulated CR spectra as well as to the H 1 emissivity measurements. 


Diffusive shock acceleration theory predicts that the particle distribution function and the 
number density of protons follow a power law in momentum dN/dp oc p~ Pi . Since flux and 
density are related by dF/dE^ in = / 3 c/ 4 - 7 t dNfdE^ in where /3 = v/c, and since dp/dE^n = 1 //3, 
we can write dF/dEkin = c/ 47 t dN/dp oc p~ p2 where Ek% n is the CR kinetic energy per nucleon. 
Therefore, as noted in Dernier (2012) and Dermer et al. (2013), if the differential flux is plotted 
versus the momentum, one expects to observe a power-law spectrum. For protons and Helium we 
adopted a spectral model similar to the one used to describe BESS fluxes (Shikaze et al. 2007): 
dF/dEkin = A/3 Pl {p/po)~ P2 where p is the momentum per nucleon and p$ = 1 GeV/c. The 
factor fd Pl allows for a possible departure from a power-law spectral shape at low energy. We 
tested several forms; this spectral shape gave the best compromise between the number of free 
parameters (A, P\, P 2 ) and goodness of the fit to the heliospheric fluxes. For electron-positrons a 
similar formula did not provide a sufficient hardening at low energies since fd ~ 1. Based on the 
measurements of AMS-01 we chose the following form that allows a harder spectrum below 1 GeV: 
dF/dEki n = A(Ek in + (Ekin + Pi)~ 0 5 )~ P2 ■ With 10 free parameters (3x3 for the Helium, proton 
and electron-positron spectral forms plus the force-field solar modulation potential) we minimized 


Table 3. Multiplication factors of Table [2] scaled by ISM abundance. 


Projectile 



Target 




H 

He 

CNO 

Mg-Si 

Fe 

Sum 

H 

1 . 

0.366 

0.016 

0.004 

0.001 

1.387 

He 

3.68 

1.363 

0.058 

0.015 

0.005 

5.12 
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Fig. 7.— Differential local 7 -ray emissivity per H 1 versus 7 -ray energy. We calculated the contribu¬ 
tion associated with hadronic collisions by folding directly the proton and Helium spectra detected 
by PAMELA and ATIC-2 with the 7 -ray production cross-section of Kamae et al. scaled by the 
multiplication factors provided in Table [3] (red dash-dotted). We calculated the bremsstrahlung 
from the electron LIS of Webber et al. scaled by a factor 0.64 (green dash-double-dotted). The 
solid line represents the sum of those contributions; it does not reproduce the LAT emissivities 
because of the CR modulation by the Sun. No free parameters were adjusted for this graph. 
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Fig. 8.— Differential local 7 -ray emissivity per H I compared to predictions based on Kamae et al. 


(2006) (black solid), Huang et al. (2007) (red dashed), and Shibata et al. (2013) (blue doted) 7 -ray 
production cross-sections. Here we apply a force-field approximation to the CR proton and Helium 
spectra with the potential of $ = 500 MV. 
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the sum of four y 2 values for four fits: the proton spectral form to a combination of heliospheric 
proton fluxes from PAMELA and ATIC-2 from 0.5 to 200 GeV, the Helium spectral form to Helium 
fluxes also from PAMELA and ATIC-2 in the same energy range, the electron-positron spectral 
form to the fluxes measured by Fermi-LAT from 7 to 200 GeV, and the corresponding emissivities 
to our Fermi-LAT 7 -ray emissivities. While the PAMELA and ATIC-2 proton and Helium energy 
ranges are relevant for 7T° production, the one of the Fermi-LAT electron-positron corresponds to 
7 -ray energies where the bremsstrahlung component is negligible. The fit to the LAT electron- 
positron measurement only fixes the asymptotic shape of the electron-positron form. We decided 
to not use any direct electron observation at lower energy to fit the electron-positron form to 
simplify the calculation of the solar modulation. The electron-positron spectrum therefore depends 
strongly on the H I emissivity as well as on the function we choose to represent its shape. For 
the evaluation of x 2 w e included systematics and statistical uncertainties in the local emissivity, as 
well as the published uncertainties for PAMELA and ATIC-2 data. The CR spectral forms were 
folded with the 7 r° cross-section of Karnae et al. (2006); we propagated an error of 20% on this 
cross-section to account for the spread between models that we observed in the previous section 
and for the energy dependence of the scaling factors that we neglected. The parameters resulting 
from the fit are given in Table [4j In addition, we find a solar modulation potential corresponding 
to the data-taking period of PAMELA considered here of (580 ± 30) MV. This value is higher than 
indirect measurements of ground-based neutron detector Usoskin et al. (2011) or of the ACE/CRIS 
spectrometer |^] which predict a modulation potential decreasing from 550 MV to 350 MV between 
July 2006 and February 2008. A possible origin for this descrepency is the under-prediction of 
the model above 1 GeV that the fit compensates for by increasing the modulation potential. As a 
consequence the experimental emissivities at low energies are over-predicted. 


In Figure [9] (bottom) we show the proton, Helium and electron-positron LIS resulting from 
the fit as well as their uncertainties. In the figure we also show the heliospheric fluxes used for 
the fit, and the electron-positron LIS. Folding those spectra with the production cross-section of 
Karnae et al. gives the 7 -ray emissivity shown as a solid line in Figure [9] (top) together with the 
la error contour. In this panel we decomposed the emissivity into a contribution from hadronic 
decay (red dash-dotted curve) and bremsstrahlung (green dash-double-dotted). We observe that a 
fit combining the heliospheric fluxes and the local emissivity under the hypothesis of the force-field 
approximation and of the LIS spectral form we chose reproduces the observations well. However 
the LAT emissivities at low energies seem systematically over-predicted while the high energy ones 
are under-predicted. The standard deviation corresponding to those variations is 14% which is 
the same order as the difference of the three production cross-sections (Figure [ 8 ]). We obtain 
an electron-positron LIS notably different from the AMS-01 observations. However, retreiving an 
electron-positron LIS from the emissivity is challenging given that the bremsstrahlung contribution 
is significantly lower than the hadronic one and that the LAT emissivities at low energy might be 
biased by unaccounted for point sources given the broad LAT PSF below 100 MeV. We note that 


%ttp://www.srl.caltech.edu/ACE/ACENews/ACENewsl55.html 
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some of the 10 fit parameters are strongly correlated; we show in Figure 10 the 2D error ellipses 
calculated at la. 


Table 4: CR spectral form parameters deduced from the fit to a combination of the LAT emissivity, 
the proton and Helium fluxes from PAMELA and ATIC-2 and from the Fermi electron-positron 
observations. 

A PI P2 

nr -2 s -1 sr -1 GeV -1 

Proton (2.38 ± 0.09) xlO 4 1.1 ± 0.6 2.838 ± 0.009 

Helium (1.07 ± 0.02) x 10 3 0.4 ± 0.3 2.777 ± 0.004 

Electron-positron (3.8 4= 0.2)xlO 2 4 4= 2 3.22 4= 0.02 


In Figure 11 we display the proton and Helium LIS versus the CR momentum per nucleon. 
Although the best fit proton and Helium LIS show a hardening at low energy, we observe that 
within the uncertainties they are compatible with a power law in momentum. In the same graph 
we plotted the LIS obtained by Shikaze et al. (2007) and by Dermer et al. (2013); the latter were 
derived from a preliminary set of the emissivities described here and with an updated version of 
the 7r° cross-section described in Kachelriefi Sz Ostapchenko (2012). 


10. Summary 

From four years of Fermi -LAT observations we have derived the 7-ray emissivity per hydrogen 
atom between 50 MeV and 50 GeV in the solar neighborhood. We have measured its spectrum 
and estimated the systematic uncertainties, and we have shown that the spectral shape of the 
emissivity is the same within uncertainties for hydrogen in the atomic, molecular, ionized or DNM 
components of the ISM. We have extracted the molecular hydrogen-to-CO conversion factor Xco = 
(0.9024=0.007) x 10 20 cm -2 (K km s -1 ) -1 and the dust-to-gas ratio Xjjust = (41.44=0.3) x 10 20 cm' 2 
mag -1 . We have shown that the WIM is not distributed like NE2001 prediction, nor does it follow 
a simple exponential scale-height. Only a small fraction of the WIM has the same morphology 
as the free-free emission. We did not observe any significant variation of the emissivity in four 
large sub-regions of the sky. We have compared the LAT emissivities with the ones predicted from 
various production cross-section models and heliospheric proton and Helium fluxes. We obtained a 
good agreement when those CR fluxes are solar-demodulated. We have fitted parametrized spectral 
forms to the CR heliospheric fluxes as well as to the LAT emissivities and deduced possible LIS 
for proton, Helium and electron-positron assuming a force-field approximation. The accuracy of 
7-ray production cross-sections was a signficant limitation in deriving the CR spectra from the 
emissivities. We have also derived the solar modulation potential corresponding to the PAMELA 
observations. 
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Fig. 9.— Proton, Helium, and electron-positron LISs resulting from the fit combining the LAT 
emissivities and heliospheric data. Top: LAT (black points) and predicted (solid line) ernis- 
sivity calculated with the Kamae et al. (2006) cross-section and proton, Helium, and electron 
fluxes represented by a spectral form. It is decomposed into an hadronic (red dash-dotted) and 
a bremsstrahlung (green dash-double-dotted) contribution. The green contour corresponds to la 
variations in the fit parameters. Bottom: The black line represents the fitted spectral form for 
proton (left), Helium (middle) and electron-positron CRs (right). The green contour corresponds 
to la variations in the CR fit parameters. When we apply the solar modulation to this contour 
we obtain the grey contour. Heliospheric fluxes for proton and Helium are provided by PAMELA 
(Adriani et al. 2011) and ATIC-2 (Panov et al. 2006| ) (above 100 GeV) observations. In the lower 
right plot corresponding to the electron flux we show data from AMS-01 (Aguilar et al. 2002) to¬ 
gether with Fermi -LAT electron-positron fluxes (above 10 GeV) but only the latter were used for 
the fit. We scaled the proton and Helium fluxes by and the electron-positron fluxes by E\ in . 
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Fig. 10.— Fit parameter error ellipses calculated at lo\ For each plot we subtracted from the 
parameter its value given in Table [4] and divided by its uncertainty. In this graph the ellipses are 
therefore centered on 0; each axis is drawn from -1.5 to 1.5. 
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Fig. 11.— Best fit proton and Helium spectral forms (black lines) versus the momentum per 
nucleon. The green contour corresponds to la uncertainty. The dashed line shows the LIS obtained 
by Shikaze et al. (2007) while the magenta shaded area corresponds to that of Dermer et al. (2013). 
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